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ABSTRACT 



A general analysis of coplanar lines on a dielectric substrate is 
presented. The field equations are written in terms of Hertzian poten- 
tials and boundary conditions are applied to generate a system of 
linear equations. These equations are then Fourier transformed and 
reduced to a form where a computationally efficient solution can be 
obtained by the method of moments. Any desired accuracy may be 
obtained by expanding the transform of the current if the problem is 
one involving strip conductors or the electric field if there are slots. 

A one term expansion produces very good results. Both wavelength 
and characteristic impedance of the transmission line structure are 
obtained. Theoretical and experimental results are presented for 
coplanar strips. 
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I. INTRODUCTION 



During the past several years a variety of coplanar microwave 
transmission lines has been developed covering a wide spectrum of 
applications. However, there still remained several possible struc- 
tures to analyze, such as, coupled- slot lines, coplanar strips, etc. 
Also, theoretical analytic methods more accurate were sought to 
reduce the errors obtained by the more classical methods, such as 
the electrostatic approximation, the method of partial images, etc. 

The present work is concerned with two parallel, coplanar strips 
over a dielectric substrate, the other side being bare, as illustrated 
in Fig. 1. The structure is uniform and infinite in both x and z direc- 
tions; it is also assumed that the substrate material is lossless and 
its relative permittivity is € 

r 

For coplanar strips (CPS) to be practical as a transmission line, 
radiation must be minimized. This is accomplished through the use 
of a high permittivity substrate, which causes the coplanar strip 
wavelength to be small compared to the free-space wavelength, and 
thereby results in the fields being closely confined to the strips with 
negligible radiation loss. 

A spectral domain transform method was suggested by Professors 
Itoh and Mittra [Refs. 1, Z] which yields an exact solution; this method 
together with the method of moments [Ref. 3], is used in this present 
work. 
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Fig. 1. Parallel coplanar strips configuration 
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II. DISPERSION CHARACTERISTICS 



A. FIELD AND BOUNDARY CONDITIONS 

Let the following electric and magnetic fields exist along the 
propagation or z-direction, as justified by the Hertzian vector poten- 
tial functions described in Appendix A, 



6* 



(Al) 

(A2) 



where IT is the propagation constant, and since it was stated that the 
substrate material was lossless, it follows that 



y = i 






(A3) 



From the field expressions in equations (Al) and (A2), all other 
components of electric and magnetic fields can be derived from 
Maxwell’s curl equations, as shown in Appendix B, as, 
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Applying boundary conditions at the interface between regions 2 
and 3, tangential electric and magnetic fields must be continuous; 
therefore, at y = 0, 
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(All) 



Similarly, at the interface between regions 1 and 2, tangential 
electric fields must be continuous and tangential magnetic fields 
discontinuous by corresponding surface current densities; therefore, 
at y = d. 
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Also, the electric fields will exist only in the dielectric part of 
the interface and can be expressed as, 
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Substituting the field expressions of equations (Al) through (A7) 

into the boundary conditions expressions of equations (A8) through 

(A17) one obtains, 
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B. SPECTRAL DOMAIN TRANSFORM 



All potential functions must satisfy the following relation: 

$ + 4c* $ “ ° 

where. 



(Bl) 



D=v ! +i 



(B2) 



One introduces the Fourier transform to the <x -domain, as 
suggested by Itch and Mittra [Refs. 1, 2], via, 
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&<«,?)-/ d * , i«U,3 i (B3) 
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Therefore, one can find the (X -domain Fourier transform of 
equation (Bl) as, 
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Furthermore, one can define propagation constants for each 



region, 
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From planar dielectric substrate theory, the effective permittivity 
is related to £ Q and G ^ as > 

(B7a) 

Therefore, analyzing the expressions for propagation constants 
in equation (B6), for region 1 or 3, 

(B8) 

where X is the dielectric's effective wavelength and is related to 
the free-space wavelength A as, 

^ < * < X . (B8a) 



Replacing this last expression into equation (B8), one can find 
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therefore, ^ ^ and ft ^ are always real quantities, independent of 
the values of 0< . 



Similarly, for region 2, 
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where, as before, the upper and lower bound values for Cx: are, 
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It is clear that 'S y will be imaginary for small values of cx and 
real for large values of . Specifically, ^ ^ will be imaginary 
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and real for 

<X> X ^- 1 ' . (B9c) 

Solving the potential functions differential equation (B5) for each 



of the three regions one obtains, as derived in Appendix C, 
field potential functions: 
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( B 1 7 ) 



c j £> (W M 

(Jj ) ( ) U 

where both, B ^ and ^ are evaluated at ^ ^ = j Y 2 

In all the following discussion, equations subscripted with lower 



(B 17a) 



case ,, a M will refer to ^ being imaginary. 

The set of equations (BIO) through (B15) can be interpreted 
physically as follows: 

Electric and magnetic fields are more concentrated in the 
dielectric substrate, region 2, and decay exponentially toward 
zero outwards into the surrounding media, regions 1 and 3. 

In the dielectric substrate itself, electric and magnetic fields 
and, hence, energy, are more concentrated in the lower values 
of OC . Although there is no direct comparison of the domains 
and X , a set of plots of the potential functions indicate 
that energy is also confined to the vicinity of the strips. 

An approximate graphical representation of the energy distribu- 
tion within the dielectric substrate for a fixed value of dimension y 
is shown in Fig. 2. 

It is convenient to obtain the Fourier transform of equations 
(A18) through (A26) as follows, 
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Fig. 2. Energy distribution within the dielectric substrate 
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where the derivative transform pair, 
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has been applied. 

Upon substitution of the field expressions of equations (BIO) 
through (B15) into the above relations, these become, 
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where, in equations (B40a) and (B41a), 

Furthermore, one can relate the constants and 

(«*) (o^r) 

in terms of |o^) and A^\ot) via equations (B32) and B33), 
obtaining, 
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A relation between the constants A^ e ^(oc ) and A^^(o^ ) in terms 
of the surface current densities J x (^£ ) and J z (o C) can be established 
via equations (B34) and (B35), obtaining, 
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Finally, applying equations (B36) and (B37), the surface current 

densities can be related to the electric fields ^ (cx) and ^§> (ex), 

x z 

obtaining, 
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where all the intermediate steps are fully derived in Appendix D. 



C. PHYSICAL PARAMETERS 

Up to this point, the problem has been stated in a quite general 
form, i. e. , there is no dependence on the actual physical configura- 
tion other than the boundary surfaces. 

In general, the surface current densities, ) and J z (^ ), can 

be expressed as a train of known basis functions, [Ref. 3], as, 

to 

Jx<«l * T> a <- (Cl) 

3* CM - Za ^ <j*i (C2) 

is * 

Substituting the above expressions into equations (B46) and (B47), 
one obtains, 

Ml («,/£> % 7a (C3) 

1 t*f U 1 ^ 

Ma(<*|S) Z M 4 Z = S*(«0 _ ( C4) 

Since the two conductors are sufficiently narrow, one may assume 
that the surface current density in the x-direction is zero, which will 
modify equations (C3) and (C4) as, 
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Furthermore, although theoretically one should perform a multi- 
term approximation to the surface current density in the z-direction, 
in this case a one term approximation is sufficient since a good 
assumption of this surface current density is a square pulse, i. e. , 
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or, substituting equations (Cl) and (C2) for the case where i = 1, 
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It is clear that one only needs to work with either one of the 

above equations; consider equation (C9), for the present case. 

As mentioned before, a good approximation to the z-direction 

surface current density is a square pulse; therefore, it is necessary 

to obtain the C* -domain Fourier transform of J (<X), shown in 

z 

Fig. 3, as, 
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Therefore, equation (C9) becomes, 
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In general, for a strip configuration, the electric fields are not 
easily defined; therefore, it is convenient to make equation (C12) 
independent of <§ z (oc). This could be done by taking an inner product 
with a function orthogonal to % (oc). Applying this concept, equation 
(C12) becomes, 
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Z - directed surface current density 
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A suitable weighting function W(oc) could be the complex conjugate 

of J (oc ), i. e. , J (- oc). 
z z 

Also, by Parseval's theorem, the right-hand side of equation 
(C13) becomes zero because of the orthogonality of the integrand. 
Therefore, equation (C13) becomes, 



It is clear that once equation (C14) is integrated, the dependence 
on the variable &c disappears, so one can state that the whole process 



teristics, namely, width of the strips W, separation between the 
strips S, thickness of the substrate D, and the dielectric f s relative 
permittivity £ . 

Once the oc -domain Fourier transform of the surface current 
density is obtained, expressions for the different constants, A through 
D, have to be redefined to account for the imaginary character intro- 
duced by J (oc). It turns out that all the magnetic potential constants, 



all the electric field constants, real. A detailed discussion of this 
procedure is presented in Appendix E. 

Using the expressions derived in Appendix E, the equations for 
the various constants are modified as, 
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D. NUMERICAL INTEGRATION 

Essentially, the main problem, as far as numerical integration 
is concerned, is the evaluation of equation (C15). Although theo- 
retically the limits of integration are - oo and + oo, truncation at 
points corresponding to sufficiently small values of the integrand is 
allowed. 

The numerical integration method used is a modified Simpson's 
rule. Consider a general function, as shown in Fig. 4, where one 
divides the axis of integration oc in N segments of width A ; the 
ordinate, corresponding to the half-width point of each segment, is 
taken as the average ordinate of each integration segment and a 
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Fig. 4. Simpson's rule of numerical integration 
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Fig. 5. Typical integrand plot 



27 



rectangle is formed with the segment^ width. It follows that the 
total area under the given curve will be approximately equal to the 
sum of the areas comprised by the N rectangles. 

It is clear that for more accurate results the width of each seg- 
ment A should be a function of the rate of variation of the function 
F(o^), i. e. , for rapid- varying functions, the width A should be 
smaller than for slow-varying functions. 

Several graphs of the integrand of the final equation (C15) were 
obtained; it is clear that for a given substrate of thickness D and 
relative permittivity €. , one can vary the frequency of operation, 

width of the strip conductors W, and separation between them S, 
obtaining different corresponding values of effective wavelengths on 
the dielectric X . Therefore, there is no single graph that could 
be used throughout the complete analysis. However, a graph corres- 
ponding to a set of parameters chosen arbitrarily will give a fairly 
good idea of how the integrand behaves and, consequently, will 
indicate the approximate limits of integration and an acceptable 
segment^ width A . A typical graph is shown in Fig. 5. 

It was chosen, according to the graphical expression of the inte- 
grand, that the limits of integration will be -1700 and +1700, since 
contributions beyond this point were less than 1% of the maximum 
contribution. Likewise, a good choice for the segment’s width 
is 0.5, which will give a total of 6800 steps of integration. 
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The rate of change of the value of the integral in equation (C15) 
with respect to the effective dielectric wavelength X was also 
investigated to determine whether it was a monotonically increasing 
or decreasing function, an oscillating function, or a non-uniform 
function. For this purpose, the integral was evaluated at several 
increasing values of X and it was observed that it behaved as a 
monotonically decreasing function, as shown in Fig. 6. 

The resulting integral in equation (C15), 

+1 Too 

J M, ~ sir? sir? | Vd d<* 

-lloo 

will reach the desired value of zero at an approximate value of 
for a given set of parameters, i. e. , frequency, width of conductors, 
separation between conductors, thickness of the substrate, and 
dielectric permittivity. 

Since ^3 is defined as 

2* 

= x # 




where X varies as 

A- c x' < X , 

v/g7 

it is clear that an iteration procedure should be used, varying 



from a starting value equal to X _, the dielectrics wavelength, 

d 



Xd - 






in appropriate steps toward X , until the value of zero is reached. 

The accuracy of the effective dielectric^ wavelength depends 
critically on the size of the above mentioned iteration steps. 
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Fig. 6. Resulting integration variation with 
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Therefore, there is a trade-off in practice of accuracy vice computer 
time or, equivalently, accuracy vice dollars. For the present work 
accuracies of the order of 0. 1% were obtained by using iteration 
steps equal to 1/1000 of the corresponding free-space wavelengths. 

The integration subroutine itself was tried for several curves 
of known equations of various rates of variations, i. e. , sinusoids, 
paraboloids, cubics, etc. In all cases variations of, at most, 0.2% 
their theoretical value were obtained. Also the limits of integration 
were increased to 2500 obtaining variations of, at most, 0. 05%. 
Finally, integration steps of 0. 4 and 0. 6 were also tried to check 
accuracy and, again, variations of only 0. 1% were observed. 

E. COMPUTER PROGRAMMING AND RESULTS 
1. Computer Program Organization 

The dispersion characteristic program, written in FORT- 
RAN IV language, accepts the following data: 

- Dielectric's permittivity 

- Width of the conductors, in millimeters 

- Thickness of the substrate, in millimeters 

- Ratio of the thickness of the substrate to free-space 
wavelength 

- Ratio of the separation between conductors to width of 
conductors . 
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Likewise, the limits of integration, the step of integration 

and the iteration step X, need to be specified. 

The program starts by calculating all the required physical 

constants, namely, H , €o > » and • All the c* - and 

p -independent parameters are then evaluated, i. e. , k^, k£, ^ , 

etc. The iteration step is initiated by assigning the lower bound 

value of X , i. e. , X , defining in this way the value of ^ . At 

this point the integration is performed by assigning the lower bound 

value of , i. e. , -1700 for the present case defining consequently 

the parameters k c , k C 2 >, k^, ft j, T 2 * an( 3 Y y Means are 

provided to test the real or imagimary character of 0 2 and two 

corresponding paths then follow. The integration, or better, the 

summation stops when c< reaches its upper bound value, i.e., +1700 

for the present case, and means are provided to test whether the 

absolute value of the resulting quantity is less than an arbitrary 

-7 

small epsilon value, i. e. , 10 for the present case. If the case is 
that the result is greater than this epsilon value, X is increased 
toward X in one iteration step X and the summation is performed 
again. The whole process ends either when a result less than epsilon 
is found, in which case the corresponding X is the answer sought, 
or when X reaches its upper bound value, i. e. , X 

The program is outlined in detail in Appendix F. 



aTjL. 
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Fig. 7. Dispersion characteristic curves for € r = 12 
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Fig. 8. Dispersion characteristic curves for € r = 16 
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Fig. 9. Dispersion characteristic curves for € r = 20 
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2. Computer results 



Effective dielectric wavelength ratio relative to free-space 
wavelength information was obtained for three different relative 
permittivities, namely, for = 12, 16, and 20. This information is 
presented in terms of ratios of the different parameters S, W, D, and 
in order to permit the use of any set of actual physical parameters. 

For each permittivity, curves were obtained for the following 
values of the ratio of separation between conductors to width of con- 
ductors, S/W = \ , % , 1, 3/2, 2, 3, and 4, corresponding to values 
of the ratio of substrate thickness to free-space wavelength D/X 
from 0.01 to 0.06. 

Dispersion curves for £ =12 are presented in Fig. 7, for 

r 

= 16 in Fig. 8, and for £ = 20 in Fig. 9. 

F. EXPERIMENTAL RESULTS 
1. Fabrication 

The coplanar strip (CPS) fabrication was made by photo- 
etching technique. Photoresist is first spread uniformly across the 
surface and allowed to dry. This film is then exposed to ultraviolet 
light through a photographic mask that permits the light to fall only 
on those areas in which the copper layer is preserved. The un- 
exposed film is next washed away from the surface, leaving a solid 
film of photoresist in the areas that were exposed. During the 
etching process the copper layer is removed in the unprotected areas 
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but is unaffected under the photoresist. Finally, the photoresist is 
removed with acetone, leaving only the dielectric surface with the two 
copper strips, as shown in Fig. 10. 

It was decided to construct the coplanar strips on a piece of 
6 11 x 3 n x 1 / 8 1 1 copper clad substrate of relative permittivity of 12. 
Also, the ratio of width of conductors to thickness of the substrate 
was chosen arbitrarily as 1 and, similarly, the ratio of separation of 
conductors to width of conductors was also chosen arbitrarily as 1. 

The coplanar strips were fed with a transition consisting of 
a 2% inch piece of 85 mil semi-rigid coaxial line, soldered as shov/n 
in Fig. 10. 

2. Experimental data 

Basically, the measuring devices used were a slotted line 
fed by a sweep frequency oscillator and monitored by a voltage stand- 
ing wave ratio meter. 

From transmission line theory, it is known that an open - 
circuited line will have voltage nulls every half- wavelength, which 
in this case will be the effective dielectric wavelength. 

Using this procedure, data was collected for a range of 
frequencies from 1 to 6 GHz and a comparison was made with the 
corresponding computer results, observing an error of 0. 1% which 
was, as a matter of fact, the accuracy that was claimed for these 
results. 
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Fig. 10. Coplanar strips fabrication 
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Table I lists the experimental results and compares them 



with the computer results. 
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TABLE I 



DISPERSION CHARACTERISTIC COMPARATIVE RESULTS 



Coplanar 


strip of characteristics 


: € =12 
r 

d = 0. 125" 
s = w = d 


Frequency ^/X from 


^/X from 


GHz 


Fig. 7 


experiment 


1. o 


0. 447 


0. 446 


2.0 


0. 435 


0. 435 


3. 0 


0. 425 


0.424 


4. 0 


0. 415 


0.414 


5. 0 


0. 406 


0. 406 


6. 0 


0. 398 


0. 399 
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III. CHARACTERISTIC IMPEDANCE 



G. THE AVERAGE POWER AS A FUNCTION OF THE HERTZIAN 
VECTOR POTENTIAL FUNCTIONS 

As stated by Collin [Ref. 4], a general expression for the time- 
average power flow in terms of the electric and magnetic fields is 

‘Rvs e z G?e[ if <§ x-U*’ a 9 ds] < G1 ) 

S 

where, 

*.a, . (G2) 

Referring back to equations (Al) through (A7), it is found that the 
equations for the electric and magnetic fields can be expressed in terms 
of Hertzian potential functions as, 
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So, substituting terms in the expression for average power, 

~ 2 [ jj ( ~ &y 












a£ 



ay a or ay M" 4 ax 



C«l* \ 

(G9) 

c / 



“(# ^ W* -ST A'# -ST ;} d * d< /J 

Taking the OC -domain Fourier transform of the average power 



expr es sion, 



*gk W) 

-QP + * «>/•><■ 4^“, ■</))(■ ■f $ <U (Kf -, j" 4 —■ 

where the Fourier transform of has been applied, i. e. , 

^ [^ ] ° "j“ * . 

From equation (G10), performing the operations indicated, 

4 ** ( 

Qt [ [( l ® (K,i{\ <& U \oc, y) -+ <& 



(G10) 



:cii) 



~cO££> 



a^ (< ) a£«> 
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It “ "/>f ^ -sf ~Jf 4 ‘“■V 1 -Sj- < G12 > 



^ +j*Y 4 ’ ^ 



. * aflfaA)* 

J */* -ay ® 

As far as the first four terms in the above expression for 
average power are concerned, they will be real. However, for the 
remaining four terms, since they correspond to cross-power terms, 
a more careful analysis has to be done. A detailed exposure to such 
analysis is given in Appendix G. 

The real part of the average power turns out to be, 
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(G13) 
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and similarly, for Y^ imaginary, 
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It is clear from the above expressions that one can integrate out 
the y-dependence of the average power, leaving only the -depend 
The integration details are shown in Appendix G, the resulting 
expressions being, 

**■ A* 

(y - y z )[ej. ( + C ie) (K) ) +f)o ( I 4 lc (k V«)| 2 )] h 2 Y t Jl 



ence. 



♦i 



+ d(»v^) [<* (c“Vo - B m m)V(|c' u wP- 1 & <m m| 1 )] 

+ (eah *V 2 |B W («)Jc“wl] ( ^r- >i ) 

+ 1>*Y Z (|& (t V«l|c“V«)-B w M|c ( V||)(-fe 2 -|^) d J 4* 

and, similarly, for Y^ imaginary. 
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Fig. 11. Z - direction surface current density 
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+p 0 (|a (M («i| 2 + 

+ ^ (~ (C lt '(«) 4 4-yUo ((C (L V4 7 + l& ((v M*)] *'V» 2r/d 

4 d (^-r; 2 )[ e *(c^- Wm)*^ (Ic w («i|M& ( V>I 1 )] 

- ZkY z " (S <e V«)|c <U (*)l -|& W (o|c (€l M)(^ 2 +lz) d } 

where in both equations (G14) and (G14a), the equality Y = Y , 

1 3 

has been applied. 



H. THE CHARACTERISTIC IMPEDANCE IN TERMS OF THE 
DISPERSION CHARACTERISTICS 



In general, one could express the average power in terms of the 
current in the coplanar strips as, 



ft* 



' T 2 ? 
- 1 Xfc 



from where, 



(HI) 



(H2) 



In the present work, it was assumed in the dispersion character- 
istic part that the surface currents through the strips were uniform; 
furthermore, the current can be expressed, as shown in Fig. 11, as 



1 4^-dt •/ (JV dx 

| d/x = VO 

where the approximation is valid since uniform surface current 
densities were assumed. 




(H3) 
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Furthermore, since it was also assumed a unit amplitude square 

pulse as the axial surface current density, i. e. , J (x) = 1, it follows 

z 

that, 



I ~ W) 



(H4) 

(H4a) 



All constants in the average power expressions of equations 
(G14) and (G14a) are completely defined by choosing appropriate 
values for the width of the strips, separation between strips, thick- 
ness of the slab, dielectric's permittivity, and frequency of operation. 

The propagation constant p is now well-defined from the value 
of the corresponding effective dielectric's wavelength X , i. e. , 






(H5) 



where X is determined from the dispersion characteristic curves of 
Figs. 7, 8, and 9. 

Therefore, the expressions to be evaluated for CPS characteristic 
impedance are, 

-f f>. ( {A <k V«)f + I I> <M ool 1 )] 

+ z ( (B Ce) M + -fyL) 0 (/£> (A V«*)| 2 '+!c ( ^C»d[ 2 )]sfnfi 2Y z d j 

+ d Y 2 i )[e 2 (c w m - b (0 m) ( I c tW oo j 2 - 1 ef a M| 2 )] 

+(<osh 2Y 2 d-l')[<s i S M (*)C | |[c (tt, C<x) J] ( -yr - V* a )j 

4 (|B ( V)|c^M -B w M|d u (K)JX4^ .*») d } d* 

and, similarly, for the case where ^ is imaginary, 
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(H6a) 



■*■ + V/)[e 2 (C fe V^ + B ( V0^(|A)|%l3<%jP)J tin 2V J 

4 a (**-**)[«, (c*>M- &“'w) +p. ((C^W-I &v M r)i 
- 2 w r* ( 3 <e, («) |c A V*t)| - 1 a A M C ce, c«))( /S?+ LI ) «t } d u 

I. NUMERICAL INTEGRATION 

The numerical integration problem, as far as the characteristic 
impedance is concerned, is the evaluation of equations (H6) and (H6a). 
Again, although the limits of integration in ex' are - oo and + oo, a 
prior knowledge of the behaviour of the integrands allows the use of 
truncated limits . 

As in the dispersion characteristic part, the numerical integra- 
tion method used is a modified Simpson's rule. Practical truncation 
limits are again -1700 and +1700 in steps of = 0.5, which will give 
6800 computing points for a given set of parameters. 

Once the physical parameters are specified, namely, the sub- 
strate's relative permittivity, the separation between conductors to 
width of the conductors ratio, and the substrate's thickness to free- 
space wavelength ratio, all the constants are well-defined except for 
the propagation constant A . However, this value can be interpolated 
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from the dispersion characteristic curves shown in part E. There- 
fore, numerically- wis e, the calculation of the characteristic imped- 
ance is simpler than that for the dispersion characteristics since no 
iteration procedure is required. 

J. COMPUTER PROGRAMMING AND RESULTS 
1. Computer program organization 

The characteristic impedance program, written in FORT- 
RAN IV language, accepts the following data: 

Dielectric’s relative permittivity 
Width of the conductors, in millimeters 
Thickness of the substrate, in millimeters 

Ratio of the thickness of the substrate to free-space 
wavelength 

Ratio of separation between conductors to width of 
conductors 

Ratio of effective dielectric’s wavelength to free-space 
wavelength 

Likewise, the limits of integration, and the step of integra- 
tion need to be specified. 

The program starts by calculating all constants and para- 
meters and directly evaluates the impedance integral; basically, the 
procedure is the same as in the dispersion characteristic part except 
for the iteration search for X since now it is completely defined. 

The program is outlined in detail in Appendix H. 
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2 . Computer results 






Characteristic impedance information was obtained for three 
different relative permittivities, namely, for C r = 12, 16, and 20. 

As in the dispersion characteristic part, this information is 
presented in terms of ratios of the different parameters, i. e. , S/W 
and D/ X , in order to permit the use of any set of actual physical 
parameters . 

For each permittivity, curves were obtained for values of 
S/W equal to %, 1, 3/2, 2, 3, and 4, corresponding to values of 
D/ X from 0. 01 to 0. 06. 

Characteristic impedance curves for £ = 12 are presented 

in Fig. 12, for £ = 16 in Fig. 13, and for Q. ^ = 20 in Fig. 14. 
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Fig. 12. Characteristic impedance curves for 12 
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Fig. 13. Characteristic impedance curves for G.- 16 
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Fig. 14. Characteristic impedance curves for 20 
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IV. APPLICATIONS 



The present method is applicable to problems based on any 
coplanar strip-type configuration and also to problems related to 
slot lines, such as coupled slots or coplanar waveguides. 

For the latter case, the configuration is, in general, as shown 
in Fig. 15, and the main difference with the coplanar strip problem 
is the knowledge of the electric fields instead of the surface current 
densities . 

Referring back to equations (B46) and (B47), these can be 
expressed as , 






(LI) 



where, 



M ( 6^1 




(L2) 



therefore, 




(L3) 



Similarly, 



JxM = 



M, (<* ( I '<,) £3 («) 






(L4) 
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and, 






(*) = JxC^) 



(L5) 



In general, for any slot-type configuration, z-directed electric 
fields can be neglected and the transverse, or x-directed electric 
fields can be approximated with square pulse-type distributions. 

Unlike the coplanar- strip problem, here it is possible to obtain 
an odd and an even type of electric field excitations across the slots. 

In either case, equations (L3) and (L5) reduce to, respectively, 



M z (w,fS) 






(*><) 



(L6) 



and 



M, 6*, fl) ^ *1 / » 

77 ■ “7 

M 



(L7) 



where i = e, o depending on the type of excitation, even or odd, as 
shown in Fig. 16. 

As far as the odd case of excitation is concerned, ^Massumes 
a form similar to J z (04) in equation (Cll), i. e. , 



<§X 0 (<*) = ~7 €,'<0 “ (s+vj) S/O § KJ 



(X 



(L8) 



Performing the necessary Fourier transform and algebraic 
reductions becomes, for the even case of excitation, 



>*e 



(k) = ~ C.OS — (s+vo) S/O ~ KJ . 






(L 9) 
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Applying Galerkin’s method and by Parseval’s theorem, in the 
same approach as for the coplanar strip problem, the final iteration 
equation for the effective dielectric’s wavelength becomes, considering 
only equation (L7), 



•foe 

;■ 



M i fa, (S) 



(*) &X. M db c = o 



(L 10) 



or. 






(Lll) 



where ) through M^( ) are defined in equations ( 122) 

through (125), respectively. 

Remaining procedure is correspondingly similar to that for the 
coplanar strip case. 
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Fig. 15. Coupled slots configuration 
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Fig. 16. Even and odd transverse electric fields 
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V, CONCLUSIONS 



The method used in this present work succeeded in providing an 
accurate set of results for the dispersion characteristic and the 
characteristic impedance of a Coplanar Strip transmission line. 

It was also discussed that an extension of this method will allow 
corresponding results for the Coupled Slots or Coplanar Waveguide 
configuration. However, for the latter case, care must be taken in 
adapting the equations to the even case of excitation due to the fact 
that the Fourier transform of the transverse electric field will now 
be real in character, instead of imaginary as in the odd case of 
excitation. 

The amount of algebraic work can be reduced substantially by 
the use of matrices, which for the present configuration will have 
dimension eight. However, the algebraic approach used here gives 
much more insight to the problem and also reduces sharply the 
required computer time. 

As far as the theoretical development is concerned, the only 
approximation used in the representation of the axial surface current 
density, i. e. , a square pulse, and the value of zero assumed for the 
transverse surface current density. Depending on the desired accu- 
racy, more complete, and indeed complex, representation of these 
surface current densities can be used. 
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In the numerical integration part of the problem, the integrating 
method itself has been proven to be sufficiently accurate. In fact, 
the only significant approximations used are the choice of the trun- 
cation limits, the integration step and the iteration step. However, 
several different sets of these choices were tried and, at least, for 
the present configuration, negligible variations were observed. 

It is clear that the two programs used, i. e. , the dispersion 
characteristic and the characteristic impedance, could be lumped 
together into one program reducing in about half the amount of compu- 
tation and, therefore, the required computer time. 
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APPENDIX A 



AUXILIARY VECTOR POTENTIAL FUNCTIONS 

A. MAGNETIC HERTZIAN VECTOR POTENTIAL FUNCTION 

Consider a homogeneous, source-free, isotropic region; hence, 
there is no charge density and 

. (l) 

From Vector Algebra, it is known that the divergence of a vector 

is zero if the vector is, in turn, the curl of another vector; there- 
fore, one can state, 

V*\7>TJ^ = O (2) 

fga-jtojuVxKh ' (3) 

For time-varying fields 

~ <!§ 

so applying equation (3), 

i? \7*TL (, 

« V ( TC h + 

I (6) 

where, 

if} * co 

Similarly, for time- varying fields, 
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(7) 



V* & = 

so, applying equation (3), 



^ 7x (~j c0 /° =-j«^o(4^ fy, +V(f>) (8) 

~ V z Hh = \? T^h + ~V<1> 

(9) 

Up to this point 0 has not been defined; arbitrarily choose 

*-**■. - (10) 

Then, equation (9) becomes, 

Also, it is known that, 

\ 7 * B ~ V* ^ ^ 

/ u^.(/!c 2 R h + v0) = o (12) 

k* V-Tq, -f ^7. V<£> » o 

Summarizing, for TM modes, 



<§ = -Jco^) V^x (C^ 

^4 - k 2 l(~h + v7^7* 

- K7x^x7Lf, 

Similarly, for TE modes, 



H ~ j«J€ Vx It-h 
%- k* lLe-+ V?^7- 7U> 

Ofc Ote 

e v7x\7x rCe. 



(15) 

(16) 
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B. TE AND TM MODES FROM VECTOR POTENTIALS 



1. TE Modes 




2 - te f' * i- 


(17) 

(18) 




where TL h satisfies: 




v 2 1^ h - f Jfcl H h - ° . 


(19) 


For TE modes, % =0; therefore, 






(20) 


from where, 




L _ h 

JCk = Ity = 0. 


(21) 


Let: 




rc h -^e**a,. 


(22) 


Then, from equation (19), 




4> h = o 


(23) 


where, 




y = v l +V . 


(24) 


Summarizing, 






(25) 


i£( =±VG ±W Ct 0 h 
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(27) 
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2 . 



TM Modes 



*= \7 xt7* 71? 

^4 \7* ^ 

For TM Modes, ^ = 0; therefore, 

z 

n?~<p e (*,y)<£ n cL^ 

and, 

^ 4> € +kc 4> e = ° 

where, 

° Y 2 4- 

Summarizing, 

*4fer™ 

lt=t rv^e*** 

•Hi = + — y ' ^ ^ St 



( 28 ) 

( 29 ) 

(30) 

(31) 

(32) 



(33) 

(34) 

(35) 



Equations (Al) and (A2) correspond to equations (33) and 
(25 ), respectively. 
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APPENDIX B 



TRANSVERSE ELECTRIC AND MAGNETIC FIELDS 



Consider Maxwell ! s point form of Ampere's law, 




"5*5 - joe 4 


(36) 


< Ma7"3? J +S y( - 3 ? '•sr) +a i(xr“"5y J 


(37) 
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from where, 
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and, due to the z-dependence of the fields, i. e. , S 


, one can 


apply, 
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(41) 


obtaining from equations (38) through (40), 
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(43) 
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Similarly, starting from Maxwell's point form of Lenz' law, 



KJ * Jg = ^ 






= -j<Op (My.G.K'+Mij CUj -tHz CLz) 



from where, 



( 45 ) 



(46) 
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(48) 

(49) 



x 



Substituting the value of (o from equation (48) into equation 
(42), one obtains, 

. ,/ j<o€ / 






^ ’ 1 r 






(50) 



Jj J_ /v ^ • ,* \ 

*ru( w~ J ^ ' • 

Similarly, substituting the values of <§ from equation (47) 
into equation (43), of from equation (43) into equation (47), and 

of ^4 from equation (42) into equation (48), one obtains, 

y 

- For Mr , 
x 
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r (1^ ’'d‘7 >**) 
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Equations (A4) through (A7) correspond to equations (54) 
through (57), respectively. 
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APPENDIX C 



DERIVATION OF FIELD EQUATIONS 



The differential equation defining the field functions is, as 



expressed in equation (B5), 




( 58 ) 



In turn, the propagation constant is expressed as, 




(59) 



As stated before, the propagation constant for air medium, 
regions 1 and 3, is always real. However, for the dielectric medium, 
the propagation constant will be imaginary for small values of and 
real for large values of <X . 

For region 1, the field function differential equation solution is 
as follows, 



Clearly, for infinite positive value of y, the field should vanish, 

or, 




(£ 2 -V) 5,6s) =o 



( 60 ) 






( 61 ) 
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Similarly, for region 3, 



3> s (*,</) - V/ §3 Ky> * o (62) 

In this case, for infinite negative values of y, the field should 
vanish, or, 

r- * 1 (63, 

sE> 5 («,y) = g* 4 /. 

Finally, for region 2, there are two solutions, corresponding to 
the real or imaginary character of the propagation constant; the solu- 
tion corresponding to the imaginary case will first be sought as 
follows , 

(«,y) + V* 5> 2 («,y) = o 

(54) 

6 - 4 jv* 



§ z (<x,L{) = Ks +K4 Crfs r z y 

For values of ck, greater than ^ V^r~ I , the propagation constant 
will be real, or 

Yzij + JCpcosti y t y ( 65 ) 

where, 






( 66 ) 

(67) 



One can define the following constants, 



69 



( 68 ) 



AM = k .2 

B M = Ks 
3, («) - K-) 
CM* 

( k) « IC^ 



then, finally, applying the functions to the electric and magnetic field 
cases, 



= A'VxfM 


(69) 


c«, < y) « 3 w, (6<) s/rfh y^cj + c <c V*) cosh y^cj 


(70) 


sin Y*.tj -f C (c, 6»<) cd5 > 2 ,ij 


(70a) 




(71) 


V» a t 


(72) 






1 1 


(73) 


5 2 (<k| (*, 1 ) = B^W smhY*cj + C <kl (K) c<jsK Y*y 


(73a) 


5 2 ^CK,y)= S>/^(«) s/V) Yiij + C^V*) ctfs Ky 


(74) 


(*></) = D (0 <*) e Y3 1 





Equations (69) through (74) correspond to equations (BIO) through 
(B15 ), respectively. 
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APPENDIX D 



MATRIX CURRENT AND FIELD EQUATIONS DERIVATION 



From equations (B32) and (B33), one could express the constants 

( K ) ( e ) (h) 

A v (k) and A ( be ) in terms of C' («) and C ( Cn) as, 



-c'Vi 






-C <u w[((«») ! [(^)“-i] ! -i- 4cup.n) s,„h r t d 



+ "A rj (fe)* 0 "' 1 r * J ]} . 

Define the following constants: 

^ { M’ KsH )*’ 0 * “A^] slnh r ‘ d 



Therefore, equations (75) and (76) become, 

a“' (k) = o, <l w (k) - kj z (*) 
a' m * »j 3 d (KI (k) 



(75) 



(76) 



(77) 

(78) 

(79) 

(80) 



(81) 
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(82) 



from where one could solve for G \ C* ) and ^ \ ) in terms of 

C \ oc ) and \ <X ) as follows. 



d (e, C«l = |"k ) 4 A (£ V*> - Kj z A {k \«)] 

fJ.Oa-OzOs L J 



i [o^Wu, 

kJl OiJ ~ Oj.0,3 J • 

Once again define the following constants: 
kJu 



9 = 



9 2 = 



9a- 



N)/ 

K)z 

KJl K)c^ ~~ n)2.^0a 

kJ i K){| — x ^ & 



Pn- ^ 



(83) 

(84) 



(85) 

( 86 ) 

(87) 

(88) 



k), f0 4 - .OijOj, 

Therefore, equations (83) and (84) can be expressed as. 



Q ce) M = P, A (£l C«) - A <M (k) 

C f V)--f> A ( V)-P y A ( V) 



(89) 

(90) 



which correspond to equations (B42) and (B43), respectively. 

Substituting equations (89) and (90) into equations (B34) and (B35), 
one obtains, 

(^ | (fer? u ) ^ ^ 

4- ^ cosh Y 2 c{ 4 kc* )" 



( 91 ) 



/(. “ Y \ r 


/ kc* y 


jlR+A 


( bi ]* r ) 




Skc3 J 


J J ' Yz 





+ c o&h Y a d | * ^>c(M 
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(92) 



- ’] p 2 

•> y % (-fer P.3 >•.«» 

♦ A " w ^ h r ***[J^ KfeN]P.W,(^J*ll 

j+c^sh y 4 d [jwyf> ~jcoe,,y 3 (hi ) 1 fj 

+ o* tY K-fei) 2 - , ]‘^]-j a>e * y ‘* }**JiW 

Define the following constants: 

«■ ■ «-? ( «&*- • j ■ * +% (fer a ) ^ ^ 

*+ 'kc! fr* cosli V t d + /\cc, 2 

®» - {( ^ K te )*- ‘1 1 ( & J* s ) “» h ^ «■ 

+ ^ costa Yz d | 

®3 = j l^h r t d [^r,? 2 . 10L [(fej) J -j] ? 2 



- «v 



<*)p, 

y* t \ccz \i 

. r z 



m) 1 *] 



^ Y*d [<oe.Y 5 (^)» 5* -«YP U -nrK^MJP,] 

-*y } 

^ • 

a «*j *«* ^ KfeN j ^ * * y I (fen 



"Opo 

- ^06 2 Y a j 



+ coshv t d[t<r? l -a;€.y3 (^£)'? l + b <r[(hi^-iJ 

}• 



(93) 



(94) 



(95) 



(96) 



- w^r, 
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Therefore, equations (91) and (92) become, 



<2, - Q* = J x (*\ 


(97) 


A, W M + Qo- 


(98) 


from where one could obtain expressions 


(e) ( h ) 

for A' ’( * ) and ; ( <x ) 


in terms of J ( cx ) and J ( &< ) as follows, 
x z 




&(*)- [q, J* (*) - Q 3 J x (e*)] 


(99) 


A (®0 - [Q-j \Jl (®<) + \Jx ( ^) ] 

Q,G?4+QzG? 3 l 2 4 J 


(100) 



which correspond to equations (B44) and (B45), respectively. 

Substituting equations (99) and (100) into equations (B36) and (B37), 



and solving for % ( ) and ( ) in terms 

z X 


of J (^ ) and J (ctf ), one 
z X 


obtains : 




— [<2, J 2 («)- Q 3 J x («ol = 

QiQn +Q 2 Q 3 ^ 


(101) 


(\j* Y [<3i M Q3 Jx(«)1 

Q,Qu+6?zQ3 L 0 L J 


(102) 


4ja)y0oY, [<5 a Jx (*) ] } c 

Define the following constants: 


(*) 


.. UTo> 

M, = 

Q 1 Q 4 +QzQa 


(103) 


<J?iG?4 -i-QzQa 


(104) 


Ms = - (cdOoY, Qi-aQ,) 

*+ QzQs ' 1 


(105) 


= (coo* Y, Qu + « V qA 

G?q. + Qi Q 5 X r . 


(106) 
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Then, finally, equations (101) and (102) become, 



Mf ( € < f (v) 4 Mz ( 

Jj( 6 <) -f M 4 (<x,f£) OxU) •“ < ox(‘<) 

which correspond to equations (B46) and (B47), respectively. 



(107) 

(108) 
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APPENDIX E 



CURRENT AND FIELD EQUATIONS RE-DEFINITION 



From equation (A3), the character of the propagation constant 
was defined as lossless, i. e. , 



*=jf . 

So, equations (77) through (80) are modified as, 



Mfc)' 






«c/3 r/«c?\2 
IA £cs/ 


- i J s-i/ifi Vi cf 


% 

It 


/ Ac? \2 f / 

\ -(ccsi [' 




dq a — 

copo r, 






+-/>•*( e 


J COS Jo Vi d j" 



Consequently, equations (85) through (88) become, 

p . ^ 

p _ • i^i 

K), k)4*t 

p dS j [^i__ 

K)i K)u + lOzfKj^j 

P w = 

k), fJ4+ | O 3 I # 



(109) 

(110) 
(HI) 
( 112 ) 

(113) 



(114) 

(115) 

(116) 
(117) 
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Similarly, equations (93) through (96) are re-defined as, 



®' = “ fcr * d 
+ ^4 006 ^ Yz d + Accf 

®*-j K&M * 4 £ (fe ) 2 *<• 

+ j^| co&b V 2 d } 

<5,= ™ b vid y | (fe) 1 |*| - <0 z Kjg r-i] IPJ } 

-«sUSJ (-Ei) 2 |P 2 l---/iRj [(e|)-i}} + °</S 

4 i*K{g)-'BI 

Finally, equations (103) through (106) become, 

kef ®i 



M,‘ + 



KA Z 



J |q*|q*-QiM 

J_ -kef (Sb 

j |<52.| < 33“ ( 5/|Q4l 

J 



M=>=j : j T («>A>Y, (cfel-v/SQ,) 

° |6?rlQ a -Q.|Q4i V ' ' 1 

V, |6?4l) 



M 4 = j 



lQz|<33 -Qi I <2*1 1 



( 118 ) 

(119) 

( 120 ) 
( 121 ) 



( 122 ) 

(123) 

(124) 

(125) 



Furthermore, since from equation (Cll) J (** ) is imaginary and 

z 

J (^ ) is zero, the expressions for the different constants A through 
x 

D need to be re-defined. 
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(126) 



Equations (B44) and (B45) become, 



P?M 



Qi 



lQi\c?3 - Qi 



|3 5 o*)I 



A-)- , — r l a. m| 

|Q 2 1Q3-G!i|QhI 

which correspond to equations (C16) and (C17), respectively. 
Similarly, equations (B42) and (B43) become, 



(127) 



(128) 



<L (e, (fc) = P>a (c, (k) 4 (P^j A, (t V«)| 

M =j [|%| | A W fa|] (129) 

which correspond, in turn, to equations (C18) and (C19), respectively. 
Finally, equations (B38) through (B41) become, 



”j 4* It KkM am ^ (£r l‘ w «lf 

C‘V] 



(130) 

(131) 

(132) 

(133) 



which correspond to equations (C20) through (C23), respectively. 
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APPENDIX F 



DISPERSION CHARACTERISTIC PROGRAM 



THIS PROGRAM SOLVES THE DISPERSION CHARACTERISTICS OF 
CGPLANAR. PARALLEL STRIPS ON A DIELECTRIC SUBSTRATE. 

PROGRAM ASSUMES FOLLOWING PARAMETERS ARE KNOWN 

- THICKNESS OF THE SLAB, DD, IN MILLIMETERS 

- WIDTH OF ThE STRIPS, WW, IN MILLIMETERS 

- DIELECTRIC'S RELATIVE P F R M I TT I V I TY , E R 

- RATIO OF SEPARATION BETWEEN CONDUCTORS TO WIDTH OF 
CONDUCTORS, SOW 

- RATIO OF SUBSTRATE'S THICKNESS TO FREE-SPACE 
WAV EL ENGT H, DRAT 10 

THE METHOD USED IN FINDING THF EFFECTIVE WAVELENGTH IN 
THE STRUCTURE, LPR IME, IS AN ITERATION METHOD, WHICH 
SOLVES FOR THE ZERO VALUE OF AN INTEGRAL IN THE ALFA- 
DOMAIN FOR PROGRESSIVELY LARGER VALUES OF EFFECTIVE 
WAVELENGTH. 

THE SUMMATION IN ALFA IS AN APPROXIMATION TO AN 
INTEGRATION IN THE ALFA-DOMAIN FROM (-INFINITY) TO 
( +INFINITY) . 

THE APPROXIMATIONS USED IN THE PRESENT PROGRAM ARE 

- THE ITERATION STEP,X 

- THE EPSILON VALUE, DELTA 

- THE INTEGRATION STEP, A 

- THE LIMITS OF THE ALFA-DOMAIN INTEGRAT I ON , ALFA 
AND B. 

IT IS CLEAR THAT ANY OF THESE VALUES MAY EE CHANGED 
ACCORDING TO THE SPECIFIC CASE OR DESIRED ACCURACY. 

THE QUANTITY EO IS THE FREE-SPACE PERMITTIVITY IN MKS 
UNITS . 

THE QUANTITY MU IS THE FREE-SPACE PERMEABILITY IN MKS 
UNI TS. 

THE QUANTITY C IS THE FREE-SPACE SPEED OF LIGHT IN 
METERS PER SECOND. 

THE QUANTITY Ml IS THE STRIP CURRENT IN AMPERES. 

THE VECTOR SOW STORES THE RATIOS OF SEPARATION OF 
CONDUCTORS TO WIDTH OF CONDUCTORS WHICH ARE DESIRED 
TO SOLVE FOR. THE NUMBER J SETS THE DIMENSION OF THIS 
VECTOR. 

PROGRAM DEVELOPED BY LT(JG) ARMANDO LUNA ECHEAND I A , 
PERUVIAN NAVY. 

SUPERVISOR - PROF. JEFFREY D. KNORR, PH. D. 

U.S. NAVAL POSTGRADUATE SCHOOL, MONTEREY, CALIFORNIA 
SEPTEMBER, 1973. 



IMPLICIT REAL*4 (M,N,K,L) 
DIMENSION SOW (7) 

DATA WW/3.17/ ,00/3 .17/, ER/ 12.0/ 
DATA A/0.5/ 

DATA B/1700./ 

DATA X/.4E-4/, DELTA/1. OE-7/ 

J =7 

P 1 = 3.14159 

E0 = l.E-9/< 36.*PI ) 

MU=4.E-7*PI 
C =3 . E8 
FACTOR= 1. E9 
W=WW*1. E— 3 
0 = 005*1 . E-3 

READ (5,100) (SOW( I ) ,1=1 ,J) 

100 FORMAT (7F5.3) 

24 READ (5,18) DRAT 10 
18 FORMAT (F6.4) 
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Do ooo ooo oooo 



F REQ=DR AT IO*C/D 

FRE=FREQ/FACTOR 

LPRI ME= C/ ( FREQ^SQRK ER ) ) 

0MEGA=2 .^PT FREO 

OE=OMEGA*EO' ER 

OM=QM EGA^MU 

K1=CMEG ASSORT (MU'- EO ) 

K 2 = K1^ SQRT(ER) 

K3=K1 
K11=K1**2 
K22=K2**2 
DO 101 I=1»J 
S=S0W( I ) 

22 BETA=2.PPI/LPRIME 
B2=BETA*”2 
KC 1 SQ=K 11-B2 
KC 2SQ = K 22" B 2 
KC3SQ=KC1 S Q 
KC4=(KC2SQ/KC3SQ) 
KC5=KC4-1 .0 
KC6=KC5** 2 
M 3 = 0. 0 
AL FA= — 1 700 . 

1 ALPHA=ALFA+A/2.0 
A2=ALPHA**2 



PHYSICAL CONFIGURATION DEPENDENCE AS DERIVED BY THE 
FOURIER TRANSFORM OF THE CURRENT DENSITIES. 

F1 = ALPH A - - (S+WI/2 .0 
F2=ALPHAv W/2.0 
S1=S IN ( FI ) 

S2 =SI N ( F2 ) 

S 3 = S 1 * 2 

S4=S2v-i- 2 

AB = ALPHA* BETA 

G1=SQRT(A2+B2-K11) 

G3 = G1 

G 1 1 =G 1 ^ * 2 

G22=A2+B2-K22 

IF ( G22 . LT .0.0) GO TO 7 

HYPERBGLIC CASE, GAMMA 2 IS REAL 

G2=SQRT (G22) 

G2D=G2-D 

SI=SINH(G20) 

C0=C0SH(G2D). 

INTERMEDIATE ALGEBRAIC STEPS 

N1=KC4*(G3/(G2''-ER ) vKC4'SI+C0) 

N2=KC4<- AB/ ( 0E--G2 ) *KC5"S I 
N 3 = AB / ( 0M*61 ) * KC 4* G3/ ( G2- FR )* KC5* SI 
N4 = l ./ (0M*G1)*( (0M*G2-<A&**2)*KC6/(GE*G2J )-SI + 
10M*G3*KC4+C0) 

N5=N1*N4+N2‘N3 
P1 = N4/N5 
P2=N2/N5 
P 3=N3/N 5 
P4=N1 / N5 

Q1=KC2SQ* ( ( G3 /G2 ) * KC4'-' P4- AB / ( CM* G 2 ) v P2 * K C 5 T S I+KC2 S Qv 
1P4* C0+KC1SQ 

G2 =KC2 SQ*" ( ( AB / ( 0M'-'G2 ) '' KC5^ : P 1+ ( G3/G2 ) :- KC4- ; P3 )T S I +P 3TC0 J 
Q 3= ( AB* (G3/G2 ) * KC4< P4-0E ; -G2’ P2- ( A8*' r 2 ) / ( CM r; ' 02 ) : ‘ KC5 
1P2 )*S I- ( ( OE/FP )=‘-G3'KC4- ; 'P2-AB-P4-ABTP4‘KC5 ) ”-CO+AB 
Q4 = ( ( AB v t' 2 ) /( CM, G2 ) - K C 5 ’• P 1 + A 3 " ( G3/G2 )~ P3« KC4+0F * G2* 

1P1 )*SI + (AS iP?+( OE/ER K-G3* KC4' P1+AB"KC5> P3 ) % C0+ ( GE / E R ) • 
2G1 

GO TO 9 

TRIGONOMETRIC CASE, GAMMA 2 IS IMAGINARY 
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non ooo 



7 G22=ABS (A2+B2-K22) 

G2=SQRT (G22) 

G2D=G2*D 
IS I=S IN (G2D) 

TC0=C0S(G2D) 

INTERMEDIATE ALGABRAI C STEPS 

N 1 = KC 41- ( G 3 / ( G 2 E R ) ? KC 4* T S I + TC 0 ) 

N2 =KC4* A B/ ( 0 E-” G2 ) *KC5 ;,; T S I 

N 3 = AB / ( CM-*G1 ) *KC4*H'<3/ (G2< ER) ’ KC5> TSI 

N4= 1. / ( QM*G1JM ( -OMr G 2- ( A3 >"*2 ) K.C 6 / ( GE **62 ) )*TSI + 

1 0 M-v G3 # K C4 ^ T CO ) 

N 5=N 1* N4+N2? N 3 
P 1=N4/ N 5 
P2=N2/N5 
P3=N3/N5 
P4=N1/N5 

01 =KC2 S 0- ( ( G3 /G2 ) *KC4* PA-AB/ ( OM* G 2 ) t P 2 : *' K C 5 ) 3 'T S I + 
1KC2SQ*P4’< TC0+KC1SQ 

Q2 = KC2S 0 ( (Ab/(GM>G2) *KC5” P 1 + ( G 3/ G2 ) * KC 4* P3 ) <r TSI + 
1P3*TC0) 

Q 3 = ( AB-< ( C- 3 /G 2 ) *K C 4 r P 4+ 0 E* G 2* P 2- ( A B* v 2 ) / ( 0 M* G2 ) ’• KC 5* 

1 P2 1 ’-TS I - ( ( OE/ EPv ) G3’-K C4 '' P 2 - AB 1 P 4- AB' P4 : - KC 5 ) ~ TCC + A B 
Q4=( ( AB*--2) /( OM*G2 )•? KC5 Pl+AB- ( G3/G2 )•' P3-- KC4-0E - G2* 

IP llvTS I + ( AB*P3+(0E/ER B-G3-KC4* P l+AB^KCE'- P 3 ) - TC0+ 

2 (CE/ERIv-Gl 

9 Q5=Q2*Q3-Q1* Q4 
M 1=KC1S 0*0 1/Q5 
M2 = 16 . ' MI/A2VS3TS4 
M3=M3+M2^A 
ALFA=Al FA+A 

IF (ALFA.GT.B) GO TO 2 
GO TO 1 

ITERATION STEP 

2 IF (M3. LT .DELTA) GO TO 19 
LPR I KE = LPR I ME + X 
GO TO 22 

19 WRITE (6,191) S, FREQ, LPRIME 
191 FORMAT ( « »,/, 3X, • LAMBDA PRIME FOR S = *,F10.5,2X, 

1 * AND FREQ = • ,E16.7,2X, ‘IS SE16.7,/) 

101 CONTINUE 
GO TO 24 
23 STOP 
END 
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APPENDIX G 



CHARACTERISTIC IMPEDANCE RELATED PARAMETERS 



It is convenient to recall the Hertzian potential functions for the 
three regions of analysis. The electric field functions are, 



(*»</) = A (e) M e -Y *^" d * 

<*,<{)= B ( %)sinh Y z cj + cosh V 2 y 

u, y) = D (e> M e Ya y 



(134) 

(135) 

(136) 



from where the corresponding y-derivatives can be obtained as, 

(137) 



3>h (i) < e \ -rify-d) 

s -Y,A (e) M e 1 



a& 2 (e ’ 

T 

■ 



= Y 2 [s (e V*) cosh V*y + &in u K<j ] 



- 3 ^ (e) r 3 < / 

3 - y 3 d m e ; 



(138) 

(139) 



As stated in the .dispersion characteristic part of the present 
work, specifically in equations ( C 1 6 ) through (C23), the magnetic 
field related constants have imaginary character. Therefore, it is 
convenient to modify the magnetic field Hertzian potential functions 
to account for this fact. 



if] *j [f& (U C«)| s.YUy + |c ( V)|c<j5h ny] 



(140) 

(141) 

(142) 
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and, again, the corresponding y-derivatives can be obtained as, 



a $, 110 

^ ' 
*f 



=j%[|B ( V)| V*y ] 



d J^L 

*f 




(143) 

(144) 

(145) 



Since an integration in the y-domain comprises three well-defined 
regions, one should separate the y-integral into three complementary 
integrals, one per each region. 

Therefore, analyzing each term in the average power expression 
of equation (G12), one obtains, 



4, SI C«, Y ) 

-\& ( %)fsinh z Y z tf + |C«<k,(W 

4 2 2> m (k) C tt, («0 si oh Y z a c«f5 ^ Y* y 

-[&>(*) ic^r 

&, (k) (*><-]) §!**(*, y) - (A <k V«^|e 2Yl( ^ ^ 

® ( z k] (K,< j) sioh 2 Viy 4 1 C w (oc)| l c<35l ^Yzij 

4 2 | tf u («)||c (U fa| 6in K Y t y cvshYzLj 



<2>y <^y 

^ 6“'wr^y^ + |c" l M| i ,i„^ny 



+ 2S C< ’(«) C (c, C*) Siob Y z ij c^sli Y*y ] 



(146) 

(147) 

(148) 

(149) 

(150) 

(151) 

(152) 

(153) 
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( 154 ) 





*< 


<32?/° 




*1 






a$f u 




*1 " 


a$i u 




d 1 





4 7l£ ( S<)!lC (W («>K'rtli 'Yt.lj o«h Y t y ] 



(155) 

(156) 



(157) 



j - x a“'m |a w (ki| 

j = y 2 [£ K, <«)|c (6 V«)|smh a ny + ib a, ^)|c w m 

+ [& (e] (t<\ |B (k V^l +C (<I («) |c (k M) Y z cj 035^1 Wy ] 
j («*«/) ^ )= v 6 ^ (< ’^)|b w («)| 

1 o*-j 
1 

j i£Vot,«f) =-Y i [|B (t V«)[c (<l f«) ossh l Y»y +B w fu)|c rt, w|smh 2 V»y 

+ (B <e, <«) {B ft V«>| + C ft, <K)|C A V«)|)&nh Y*y c<?sh Y z ij j 

j &?'(«. l) = -fa D W («c)|D a VjI 



j « - r, a“Vk) | A w M I e' 2Y ' d) 

g <3y 1 

• ^ ■ Y Z [B) (c) H|c (0 ^|c^sh z Y*y +|B>V)|c ft V«) s»V>tf 11 y 

+ (S (< V«)|B <M <»o| + C (C '(*) |c (k Voc)|) S/'nh Y*y cosh >iy ] 
j ^%^ / c f ) = Y 3 ^Vo<) j/W*? 



J dy 



(158) 

(159) 

( 160 ) 
(161) 

(162) 

(163) 

(164) 

(165) 

(166) 
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( 167 ) 



j = Y, A (c '(^l^ U C^)|e' ZV ' <l f' dl 

v 1 £tj 1 

j -v 2 rB^V«) 1 c^Vk)| c^y^ 4 1 ^nh z \o, 

1 ' (168) 

+ (B^lB^I + C^lc^Vool) &nh Y z j ojsIi r z cj J 

j t.%) = - v, B^)| e 2 ^ . « (169 > 

Therefore, the expression for average power in equation (G12) 



becomes, 

( M ) . J- jj Y, 7 )(e, I tf'MPle- 2 ' 5 '’-" 

+ <*>/& ( *’-Y a *)(e. Ic«(«)| ! ) e^1 

iinh'ritj (€*|& <fl Ml*+ytf.|B tfc W) -«t^V, l (€ a |c W, fK)| 4 4 / U B |C ft VK)|*) 

+0^2? Y t (|B (M («)|d (e, ^-B tt, Mlc (M C«)n + Y 2 (B £4 VK)|C (t V«)|-|B%4c <£ V)) J 70) 
Y 2 ij [x^o (e 2 |c w M 2 +/^ |c w Go| l ) - ftp Y* 2 («* | £> £< M 2 4 jj. | B (M M? ) 

+ «!*? Yz {?> m M -!2 > (U (u) lc (e, CK)) + c< Y* (f3 (k V«)|d w 6K)-^l|c (l M)] 

ainh lYzij [wJL (€* B (€, Ck) £ (c V*) 4^Po } dcxct^ 



and, for the case where is imaginary, the following relations hold, 



=jV 



O.M 



smh V 2 1| a j ^ 

^rsl^ Y 2 1| t ctfs Yj." t j 
iinh^Yicjs-sin'Yj 
Y t t| = CcJs 2 Y 2 " <j 
>inh 2Y z ij =j£in lYi y 
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Therefore, equation (170) becomes, 



->rcof> (^~y^)(^ |d (€ V«)| 2 4/jo |d {M («)P) e 2V3< / 

- sin 2 Y z ' <j (e 2 ( + fJo | B w («)| j ) 4 00 (i>r z l (c 2 jC^l^jcV) 

-•^S? v (|B> cU («)|c w f«) - & <e, C«) 1 c (U mI) 

-«<^ v (& (, V«) |c (U c»o| -|£> (k| ^)|c w w) ] 

♦ r c *^ [« J ^o (ejc^jv^lc^pj+o)^ y 8 ' 2 (€e|fe (c W 2 4/^jB%)p) 

-«(£ Yj" ( S ( *VcO !c A M - i B a Vv)|c <cl C«l) 

- e* kf Y 2 " 0& (t V«)[ C (e) M-b ( %) Jc CL Vk)I) J | dudij t 



(170a) 



These two expressions correspond to equations (G13) and (G13a), 
respectively. 

The y-dependence of the average power disappears through 
corresponding integrations on regions 1, 2, and 3, as indicated 
below: 

- Region 1, 




(171) 



- Region 2, 




(172) 



(173) 



(174) 



(175) 



( 176 ) 
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( 177 ) 



[ sin 2 Y i *y ay =-L (1-005 . 

- Region 3, 



(178) 



KVss. 

-*o 

Substituting equations (171) through (178) into equations (170) and 



(170a), one obtains equations (G14) and (G14a). 
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APPENDIX H 



CHARACTERISTIC IMPEDANCE PROGRAM 



C THIS PROGRAM SOLVES THE CHARACTERISTIC IMPECANCE CF 

C CCPLANAR PARALLEL STRIPS ON A DIELECTRIC SL3STRATE. 

C 

C PROGRAM ASSUMES FOLLOWING PARAMETERS ARE KNCWN 

C - THICKNESS OF THE SLAb,DD, IN MILLIMETERS 

C - WIDTH OF THE STRIPS, WW, IN MILLIMETERS 

C - DIELECTRIC'S RELATIVE PERMI TT I VI TY , E R 

C - RATIO OF SEPARATION BETWEEN CONDUCTORS TO WIDTH OF 

C CONDUCTORS, SOW 

C - RATIO OF SUBSTRATE'S THICKNESS TO FREE-SPACE 

C WAVELENGTH, DRAT 10 

C - RATIO OF EFFECTIVE DIELECTRIC'S WAVELENGTH TO FREE 

C SPACE WAVELENGTH, LRATIG 

C 

C THE METHOD USED IN FINDING THE CHARACTERISTIC 

C IMPEDANCE IS A DIRECT SUMMATION CF INFINITESIMAL 

C VALUES BY A MODIFIED SIMPSON'S RULE. 

C 

C THE SUMMATION IN ALFA IS AN APPROXIMATION TO AN 

C INTEGRATION IN THE ALF A-D0MA1 N FRCM (-INFINITY) TO 

C ( + INF IN IT Y ) . 

C 

C THE APPROXIMATIONS USED IN THE PRESENT PROGRAM ARE 

C - THE INTEGRATION STEP, A 

C - THE LIMITS OF THE ALFA-DOMAIN INT EGR AT I ON , AL FA 

C AND B. 

C IT IS CLEAR THAT ANY OF THESE VALUES MAY BE CHANGED 

C ACCORDING TO THE SPECIFIC CASE OR DESIREC ACCURACY. 

C 

C THE QUANTITY EO IS THE FREE-SPACE PERMITTIVITY IN MKS 

C UNITS. 

C THE QUANTITY MU IS THE FREE-SPACE PERMEABILITY IN MKS 

C UNITS. 

C THE QUANTITY C IS THE FREE-SPACE SPEED CF LIGHT IN 

C METERS PER SECOND. 

C THE QUANTITY Ml IS THE STRIP CURRENT IN AMPERES. 

C 

C THE NUMBER JJ SETS THE NUMBER OF DRATIO VALUES. 

C THE NUMBER OF IL SETS THE NUMBER OF SOW VALUES. 

C 

C PROGRAM DEVELOPED BY LT(JG) ARMANDO LUNA ECHEANDIA, 

C PERUVIAN NAVY. 

C SUPERVISOR - PROF. JEFFREY B. KNCRR , PH. D . 

C U.S. NAVAL POSTGRADUATE SCHOOL , MCNTE RcY , CAL I FCRNI A 

C SEPTEMBER, 1973. 

C 

C 

IMPLICIT REAL*4 (M,N,K,L) 

DATA WW/3.17/ ,DD/3.17/, ER/20./ 

DATA A/0.5/ 

DATA B/ 1700./ 

I L=7 
J J = 6 

PI=3. 14159 
F PI =4. vp I 
EC=l.E-S/( 36.*PI ) 

MU =4. E-7 X P 1 
C =3 .E 8 
FACT GR= 1 . E9 
W=WWL 1. E-3 
D=DD*1. E-3 
M1=W**2 
DC 11 I =1 , IL 
READ (5,101) SOW 

101 FORMAT (F5.3) 

S = sew?- w 

DO 13 J=1,JJ 

READ (5,102) DRAT IO,LRATIO 

102 FORMAT (F6.4,F8.6) 
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OOO ooo oooo 



FREQ=DR AT IO^C/D 
FRE=FREG/ FACTOR 
L AMBDA=C/FREQ 
L PR IME = LRAT 1 U*L AM BDA 
0NEGA=2.-PI- FREQ 
OE=OMEGA*EO- ER 
CN = CMEG A-MU 
K1=CMEGA*SQRT(MU*E0) 

K2=K1^SQRT ( ER) 

K3 = K1 

K 11 =K1* *2 

K22=K2-' ! *2 

22 BETA=2 .*PI/LPRIME 
B2 = BETA 
KC1SQ=K 11-B2 
KC2SQ=K22-B2 
KC3SQ=KC1SQ 
KC4=(KC?SQ/KC3SQ) 

KC5=KC4-1 .0 
KC6=KC5**2 
M 3 = 0 . 0 
A IF A=— 1 700. 

1 ALPHA=ALFA+A/2. 0 
A2 = ALPHAS *2 

PHYSICAL CONFIGURATION DEPENDENCE AS DERIVED 3Y THE 
FCURIER TRANSFORM OF THE CURRENT DENSITIES. 

Fl=ALPHA*(S+W)/2.0 
F2 = ALPH A-> W/2.0 
S1=SIN( FI) 

S2=SIN( F2) 

C J= (4./ ALPHA)-Sl*S2 
AB = ALPHA-' BETA 
G1=SQRT (A2+B2-KI1 ) 

G3 = GI 

ALGI=A2 /Gl-Gl 
G 1 1 = G 1* ■ 2 
G22=A2+E2-K22 
IF (G22.LT.0.0) GO TO 7 

HYPERBOLIC CASE, GAMMA 2 IS REAL 

G2 = SQRT ( G2 2 ) 

G2D=G2-D 
G 220 = 2. -’G2D 
ALG2=A2/G2-G2 
ALAG2 =A 2+G22- 
SI=SINH(G2D) 

CC=CQSH (G2D ) 

- S I 2 = S INH ( G22D ) 

FCO=COSh( G22D )- 1 . 

INTERMEDIATE ALGEERAIC STEPS 

N1=KC4* ( G3 / ( G2-> ER )-*KC4 -SI+CO) 

N2=KC4v AB/(0E*G2)^KC5^SI 

N3 = AB/ (0M*G1 MKC^-*G3/(G2 -ER ) ' KC5>' SI 

N4 = 1 . / ( CM* G I ) *■ ( ( GM’*G2- ( AB v<> 2 ) • K C.6 / ( 0 E- G2 ) ) > S I +0M'* GO-KC 
14*C0) 

N 5 = N 1 * N4+N 2- N 3 
P 1 = N4 / N 5 
P 2=N2/N5 
P3 = N3/N5 
P4=N1/N5 

Q 1=KC 2SQ X ( I G3/G2 ) *KC 4- ! 'P 4- AB /( OM f-G 2 )“ PZr- KC 5 ) >S I + KC2 SQT-P 
14* CO+KC 1 S Q 

Q2=KC2SQ-« ( ( AB /( 0M’ k G2 ) '-KC5 PI (G3 / G2 ) ‘ KC4 v F3 )*SI+P3*CC) 
Q 3= ( AB^- ( G3 / G2 ) •>■ KC 4 < P4-0E ; G 2*P 2- { A B** 2 ) / ( OM- G2 ) v KC 5* P2 ) 
1* S I — ( ( 0E/ ER)~ G3 EKC4-. F2-AB*’P4-AB ■ i! P4*KC5 ) ^CO + AB 
Q4=( ( AB**2) /(0MvG2) <KC5 P 1 + AB ■> l G3 /G2 ) P3 KC4 + GE 5 *G2- ; PI ) 
l*Sl + ( AB*P3 + (0E/Ek)*G3'*KC4‘ P L + AE> ; KC5-'P3) »-CC + ( OE /ER ) - G 1 
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Q5=Q2*Q3-Q1-'Q4 
AE=Q1/Q5*CJ 
AH=Q2/Q5 "-CJ 
CE=P1*AE+P2*AH 
CI-=P3*AE-P4^AH 

BE=l./(CE*G2) i ((CE/ER)i‘G3 --KC4' C E + A B-^K C5 > C E ) 

B H = 1 . / { OM *G 2 ) * ( A B 1- KC 5 * C E + 0 M-- G 3* KC 4*C H ) 

DE=KC4v CE 

DH = KC4" CH 

AE2=AE**2 

Ah2=AH**2 

BE2=BE«=*2 

BI-2=BH**2 

CE2=CE**2 

CH2=CH**2 

DE2=DE? * 2 

DP2=DH**2 

T1 = (AE2 + DE2)- ? E0 

T2=(AH2+DH2)>MU 

T3=(BE2+CE2 ) - EG* ER 

T4= ( 8H2+CH2 ) > MU 

T5=(CE2-BE2HE0 i ER 

T6= (CH2-BH2 ) >-MU 

T7=BE^CH— BH^CE 

T8=BE*C E^EO-i ER+MU*BH^CH 

T 9 = T1+T 2 

T10=T3+T4 

T 1 1 =T 5 + T 6 

P=CMEGA*BETA: (ALGl i T9+0.5«ALG2-S $ I2^T10 + 0 A T 1 1SALAG2 + FCO 
1*T8*ALG2) + 2.*D* ALPHA*G2?T7v (B2-K22) 

GC TO 9 

TRIGONOMETRIC CASE, GAMMA 2 IS IMAGINARY 

7 G22=ABS(A2+62-K22 ) 

G2=SQRT ( G 2 2 ) 

A G2=A2/ C-2+G2 
A L0G2 = A 2-G22 
G2C=G2v-D 
G 22 D=2 . ; ‘ G2 D 
T SI =SI N (G2D) 

T CO=COS ( G20 ) 

T S I2 = S I M G22D ) 

TC02=1.-C0S(G22D) 

INTERMEDIATE ALGA ERA 1 C STEPS 

N 1 = KC4T ( G 3 / ( G2-ER J < KC4 ' TS I + TC.0 ) 

N2 = KC4vAB/ (OE* G2 )*KC5*TSI 

N3=AB/( GM*G1) VKC4T G3/(G2> ER ) KC 5- TSI 

N4= I e / ( Ch’f-Gl )^( ( -GM+G2-(AtJ**2 )-*KC6/( OE’-G2 } J T S I C M-* G 3* 
1KC4«TC0J 
N5=N1*N4+N2> N3 
P1=N4/N5 
P2=N2/N5 
P3=N3/N5 
P4=N1 / N5 

Q 1=KC 2SQ" ( (G3/G2) *KC4<- P4-AB/ (GM-*-G2 )■■ P2* KC 5 ) “ TS I + K C2 S Q». 
1P4STC0+KC1SQ 

C2 = KC2SCT ( ( 4B/ ( OMT-G2 ) * KC5 ? P1+IG3/ G2) KC4*F3 J-TSI+P3-TC 

10) 

03 = ( AB< ( G3 / G 2 ) v-K C 4 P 4 +0 E *. G2vP 2- ( AB* W 2 ) / ( OM * G2 ) ^KC 5* P2) 
l*TSI-{ ( CF/ ER)'- G3^KC4 > P2-A6' P4 -A B- P4-KC5 )?TCC+AB 
Q4= ( ( AB^ - 2) /(OMrGZJ'KC 5 -:PH-Ad. ( C 3 /G 2 ) < P3*- KC4-CE*G2«- PI ) 
1*TS I+(AB=«-P3 + (0E/EP. HG3*KC4-.-Pi*A&:- KC5*P2 ) v TC C+ ( OE / ER ) =v-G 

11 

Q5=Q2*Q3-Q1=- Q-4 
AE = C1/Q5 "CJ 
Ah = Q2/Q5‘ C J 
CE=P1*AE+P2~AH 
Ch=P3*A E-P4 1 AH 

8E = l./( 0E*G2)4< ( CE / ER ) *- G3^KC4^C E + AE‘; KC5 ^ C h ) 
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BH=1./ (0M*G2 )*(AB’-KC5 * CE+OM* G3 + KC4 *CH) 

DE=KC4*CE 

DH=KC4 ;t CH 

A E2 - A E+ +2 

Ah2=AH**2 

BE2=BE**2 

BF2=BH+*2 

Ct2=CE**2 

CF2=CH+*2 

DE2=DE+*2 

0F2=DH**2 

Tl= (AE2+DE2) -EO 

T2 = ( AH2+DH2 ) * MU 

T3=(8S2+CE2)*EG*ER 

T4=(BH2+CH2)*MU 

T5 = (CE2-BE2) 4: E0^ER 

T 6=( CH2-BH2 } + MU 

T7 = BE*CH-BK' CE 

T 8 =BE* C E* E 0* E R+ M U* B H 1 '- C H 

T 9=T1+T 2 

T 10=T3 + T4 

711=15+ T6 

P=GME G A *B E TA* ( A L G 1 T9 + 0 . 5 *-A G 2 *T 1 0* T S I 2+ D-> A L C G 2 t T 1 1 ) -2 . 
lO^C" ALPHA* G2*T7 ~ ( B2+K22 ) 

9 M2=A6S( P/FPI ) 

M3=M3+M 2~A 
A LFA=AL FA + A 

IF (ALFA.GT.B) GO TO 2 
GO TO 1 
2 M4 = . V 3/M1 

WRITE (6,19) ER,S0W,DRATJ0,FRE,M4 
19 FORMAT (• ' , 3Xt 1 FOR ER = • , F5 . 2 , 2X ♦ • S/W = ',F5<2,2X,'D 

1/L = • ,F6.4,2X,' FREQ = ' , F5 .2 , 1 X , • GHZ ' , 2X , ' Z IS »,E14. 
17, IX, 'OHMS ' , //) 

13 CONTINUE 
11 CONTINUE 
23 STGP 
END 
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